Giant bubble pinch— off 
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Self-similarity has been the paradigmatic picture for the pinch-off of a drop. Ifere we will show 
through high-speed imaging and boundary integral simulations that the inverse problem, the pinch- 
off of an air bubble in water, is not self-similar in a strict sense: A disk is quickly pulled through 
a water surface, leading to a giant, cylindrical void which after collapse creates an upward and 
a downward jet. Only in the limiting case of large Froude number the neck radius h scales as 
h{— log /i)^^'' cx r^^'^, the purely inertial scaling. For any finite Froude number the collapse is slower, 
and a second length-scale, the curvature of the void, comes into play. Both length-scales are found 
to exhibit power-law scaling in time, but with different exponents depending on the Froude number, 
signaling the non-universality of the bubble pinch-off. 

PACS numbers: 47.55.df, 47.55. db,47.20.Ma 



The pinch-ofF of a liquid drop is a prime example of 
a hydrodynamic singularity and has been studied exten- 
sively in recent years H S S IE IE • It has become 
paradigmatic for self-similar behavior in fluid dynamics: 
After appropriate rescaling, the shapes of the pinching 
neck at different times can be superimposed onto a single 
shape j|4y5i, . With the exception of some pioneering 
work [8|, |9j , the inverse problem of the collapse of a gas- 
filled neck surrounded by a liquid has not attracted much 
attention until very recently, with the analysis of the 
pinch-off of a bubble rising from a needle and the break- 
up of a gas bubble in a straining flow 0, 0, 0, Q| . 
The time-evolution of these collapsing gas-filled necks is 
found to follow a power law. If the dynamics near the 
singularity are solely governed by liquid inertia, then the 
radius of the neck h expressed in the time r remaining 
until collapse scales as h cx r^/^ [EIE0jO|i or, with a 
logarithmic correction, as log ft,)'^/'' cx r^/^ De- 
viations from this exponent of 1/2 are reported to occur 
only due to the inclusion of other effects. The collapse 
may be slowed down by viscosity {h (x t UlL Ii?|) or 
surface tension {h oc t ''^ or accelerated by the iner- 
tia of the gas flowing inside the neck, leading to /i c>c r^^'^ 

Hi. 

In this paper we focus on another example of this "in- 
verse pinch-off" , namely the violent collapse of the void 
created at a fluid surface by the impact of an object. 
Here we find exponents which deviate substantially from 
1/2, even though the dynamics are shown to be purely 
governed by liquid inertia, without significant contribu- 
tions from the effects mentioned above. The self-similar 
behavior h{— loghY^^ cx r^/^ appears to hold only in the 
asymptotic regime of very high impact velocities. 

In our experiment, a linear motor is used to drag metal 
disks with radii hdisk between 10 and 40 mm through an 
air/water interface with a well-controlled, constant veloc- 
ity V between 0.5 and 3 m/s (see Fig.QJi). A rod running 
through a seal in the bottom of a large water tank con- 



nects the disk with the linear motor. This arrangement 
generates giant voids in a very controlled fashion. The 
advantage of this setup is that the velocity is a control 
parameter and not the response of the object to the fluid 
forces upon impact. Secondly, due to the large scale of 
the experiment, viscosity and surface tension play a neg- 
ligible role 13 . Therefore the only important dimension- 
less parameter is the Froude number Fr — V'^ / (hdiskg), 
the ratio of kinetic to gravitational energy, which ranges 
from 0.6 to 90. The large scale of the experiment is also 
advantageous for the observation of details during the 
impact and collapse process, which is imaged with dig- 
ital high-speed cameras with frame rates up to 100, 000 
frames per second. 

A typical series of events is seen in Fig. ^-d. The im- 
pact of the disk creates an axisymmetric void which flrst 
expands until the hydrostatic pressure drives the walls 
inward. The inward moving walls collide and cause a 
pinch-off at some depth below the undisturbed free sur- 
face. The energy focusing of this violent collapse cre- 
ates a strong pressure spike on the axis of symmetry 
which releases itself in a downward and an upward jet 
[islllil . The latter reaches heights exceeding 1.5 m for 
the higher impact speeds in this experiment. It is this 
dominating role of inertia that makes our system differ- 
ent from other pinch-off processes in the literature. At 
higher recording speeds the pinch-off can be investigated 
in more detail as in Fig.^-h. There is a clear loss of both 
azimuthal and axial symmetry in Figs.^ and^', which 
can be attributed to a combination of the same conver- 
gence effect that causes an instability in a collapsing bub- 
ble miiiiii, and a Kclvin-Helmholtz instability due 

to the rapid air flow in the neck. The latter increases 
with increasing Froude number and limits the range of 
our experiments. Another factor which limits the Froude 
number range is the so-called surface seal, in which the 
void closes at the water surface as the crown-like splash 
is entrained by the air flowing into the expanding void 
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FIG. 1: Snapshots of the formation and coUapse of a surface void for the plunger experiment: A hnear motor puUs down a 
disk of radius hdisk ~ 30 mm through the water surface at a constant velocity V = 1.0 m/s (Fr = 3.4). (a-d) The collapse of 
the void imaged at a 1000 frames per second. The white lines (overlay) are the void profiles obtained from boundary integral 
simulations with the same initial condition, and without the use of any free parameter, (e-f) Details of the collapse imaged at 
12800 frames per second, (g-h) Details of the collapse imaged at 48000 frames per second. Note that the field of view decreases 
with increasing frame rate. In the very last stages of the collapse (f-g) there is a Kevin-Helmholtz instability that complicates 
the determination of the time of collapse. Immediately after the collapse air is entrapped, both in the form of a large bubble 
above disk (d) and as a cloud of microbubbles at the collapse point (h). 



|20ll21| . This process, which occurs at large Froude num- 
bers, changes the pinch-ofF considerably since in this case 
the gas pressure inside the void differs appreciably from 
that of the ambient air. 

In view of these experimental limitations, we per- 
formed numerical simulations using a boundary integral 
method based on potential theory without ambient gas. 
There is an excellent agreement between the numerical 
calculations and the experiments, as seen in Fig. ^-c. 
Here, the numerical void profiles (the solid white lines) 
coincide very well with the experimental profiles in the 
pinch-off region without the use of any adjustable param- 
eter, either in space or in time. 

To further quantify the pinch-off process, we now turn 
to the time evolution of the neck radius ^.(t), measured 
at the depth at which the void eventually closes. Because 
both length and time scales become very small close to 
collapse, it is not feasible to experimentally observe the 
collapse with only one high-speed camera recording |25l | . 
Due to the reproducibility of the experiment, we over- 
came this difficulty by matching several data sets im- 
aged at different frame rates, increasingly magnifying the 
region around the pinch-off. Figure |2K contains a dou- 
bly logarithmic plot of h{T) (compensated with t^^^) for 
both the high-speed imaged experiments and the numer- 



ical calculations, again showing excellent agreement for 
different Froude numbers. In this graph, a straight line 
corresponds to the power law behavior h = /iqt"'' . The 
exponent ah is plotted as a function of Fr in Fig. 
Clearly, there are large deviations from the suggested 
behavior ah = 1/2. Can these be explained by a log- 
arithmic correction as proposed in 

Let us first establish the origin of this logarithmic cor- 
rection in our system. Near the neck, the fiow induced 
by the collapsing void looks very much like that of a col- 
lapsing cylinder, while it must look like that of a sink, 
plus its image in the free surface (i.e., a dipole) in the 
far region. In the language of singular perturbations, the 
former would be the inner region and the latter the outer 
region; a complete descriptions would require the match- 
ing of these two regions. If we disregard the outer region, 
we can use a two-dimensional version of the Rayleigh- 
equation, which describes the collapse of an infinite cylin- 
drical cavity under uniform pressure 0, ilH |22| 
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(1) 



The pressure difference driving the collapse has been 
equated to pgZ ^ where Z denotes the depth below the 
fluid surface, which implies that the system is composed 
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FIG. 2: (a) The radius of the void at the depth of closure h, compensated with r^^^, as a function of the time r remaining 
until collapse in a doubly logarithmic plot, for three different values of the Froude number Fr. Experiment (blue symbols) and 
numerical simulations (red symbols) are seen to agree very well for Fr = 3.4 and 10.1. For Fr = 163 only numerical data are 
presented, because for this Froude number experiments are hindered by the surface seal (see text). The error bars, indicating 
the error in the experimental data, are usually small, but occasionally become very large for frames very close to the collapse 
time, (b) Doubly logarithmic plot of the radius of curvature of the void profile R compensated with r^''^ as a function of r for 
the numerical simulations of (a). Both h and R are well described by power laws for up to four orders of magnitude in r. (c) 
Power-law exponents ah for the radius of the void at closure depth hij), ay for the radius of the void including the logarithmic 
correction 2/(t) = h{—\ogh)^^'^, and an for the radius of curvature of the void at closure depth R, all as a function of the 
Froude number. 




of non-interacting horizontal layers of fluid, with a neg- 
ligible vertical velocity component |2^. Although the 
quantity ft-oo must in principle be determined by the 
matching process alluded to before, it is expected to be 
of the order of a typical length scale of the process, such 
as the cavity depth. Thus, strictly speaking, h^a is a 
function of time and of the Froude number. However, 
near pinch-off, the time scale for the neck motion is much 
faster than that for the evolution of the other parts of the 
cavity so that hoa may be considered only a function of 
Fr. After an initial expansion of the void, the collapse 
starts from rest at a maximal radius h^^ax (of the order 
of hdisk)- Using this as an initial condition, and treating 
hoo as a constant, the energy integral of Eq. can be 
readily found: 



asymptotic regimes 




1 



\og{h/hoo) 



(2) 



where we have introduced the non-dimensional variables 

h EE h/hmax, hoo = hoa/hmax-, and f = Ty/gZjK^^. 

Close to pinch-off, h? <^ 1, such that — 1 w 
With this approximation, we can integrate Eq. ^ once 
more to arrive at 



2\og{hoo/h) 



2 f = hyiog{hoo/h) + ^^hl^ erfc 

(3) 

For small h the term with the complementary error func- 
tion is always small compared to the first one and their 
ratio vanishes for h [0. Neglecting this term we find two 



(log(/ioo)) = \/2fi/2 for /i> l/Zioo , and(4) 
/i(-log(/i))^^^ = %/2fi/2 for/i<l//ioo. (5) 



From Eqs. Q and (0) we conclude that the scaling de- 
pends crucially on the value of hoc'- Initially, for the in- 
termediate regime h 3> h'^^^/hoc, we expect to find a 
strict power law h cx t^/^, since log(/ioo) is constant. For 
times closer to the pinch-off, when h <^ h'^ax/^oo, loga- 
rithmic corrections play a role, and the power law should 
be modified into Eq. ||SJ). 

As hoo ~ hmax in our experiments, the latter inequal- 
ity can be read as ft. <C h^ax , which is satisfied in most of 
the region where /i(t) asymptotically behaves as a power 
law (cf. Fig. |2t). We conclude that in our system the 
logarithmic correction cannot be neglected. If we plot 
the quantity y = /i(— logft)^/^ vs. time, we again ob- 
serve a power law y cx r"" , but with a slightly different 
exponent ay than the one found for h. In Fig. [21; we com- 
pare ah and ay as functions of the Froude number. As 
discussed before, an is substantially larger than 1/2, but 
even if the logarithmic term is included we continue to 
find a slower collapse for low Froude numbers. Although 
the logarithmic correction does bring the result closer to 
the suggested value 1/2, it cannot account for all of the 
observed deviations. 

Clearly, the observed anomalous power law behavior of 
the neck radius must reflect itself in the time-evolution of 
the free-surface profiles of the collapsing void. If the pro- 
cess were self-similar, the free-surface profiles at different 
times r would superpose when scaled by any character- 
istic length, e.g., the neck radius h. Actually, it is found 
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that the depth of minimum radius increases somewhat 
as the collapse progresses and it is therefore necessary 
to translate the profiles in the vertical direction so as to 
match the position of the minimum radius point before 
attempting this operation. Even if this is done, how- 
ever, the results fail to collapse onto a single shape. The 
rescaled profiles become more and more elongated as the 
pinch-off is approached which proves that the collapsing 
void is not self-similar in a strict sense. 

The free-surface shapes near the minimum point 
should thus not only be characterized by ft.(r), but also by 
a second length-scale, the radius of curvature R{t) in the 
vertical plane (see Fig.^). The spatial resolution of the 
high-speed camera images limits the accuracy with which 
R can be extracted from the experimental observations, 
but this quantity is easily determined from the numerical 
calculations (see Fig. IJt). When the radial dimensions 
{H, cf. Fig. ^) are scaled by h and vertical ones [Z) 
by VhR, the profiles do collapse, which may only sig- 
nal that their shape is very close to parabolic [23|. The 
time-evolution of this radius of curvature is also found 
to follow a power law, R = Rqt""^, the exponent an of 
which increases with the Froude number as can be seen 
in Fig. 13; ^S- 

The essence of the time-evolution of the void profile 
and the departure from self-similarity in the strict sense 
is captured in the aspect ratio h/ R of the collapsing void, 
h{T)/R{T) = {ho/Ro)T^°"'~°"'\ in which the prefactor 
ho/Ro and the exponent {ah — an) both are found to 
depend on the Froude number. It is seen in Fig. Et 
that ah — an > for any finite Froude number, caus- 
ing the ratio h(T)/R{T) to vanish in the limit r ^ 0. 
This means that in this limit R{t) becomes large with re- 
spect to the neck radius, elongating the profiles more and 
more towards the cylindrical shape close to the pinch-off, 
thereby justifying the assumptions made in the deriva- 
tion of Eq. ^ in the limit r ^ 0. A numerical fit gives 
{ah — Qffl) oc Yr~^'^'^, which indicates that h and R have 
the same time dependence as Fr ^ oo and, therefore, 
that self-similarity is recovered in this limit. 

A second numerical fit shows that Hq/Rq oc Fr~°'^°, 
which tends to zero as Fr cxi. This feature expresses 
the experimental observation that the initial elongation 
of the neck is larger for large Froude number, which effec- 
tively increases the time-interval for which the assump- 
tion of pure radial flow is valid [cf . Eq. |^] . 

In conclusion, our experiments on the collapse of a gi- 
ant surface void are in excellent agreement with bound- 
ary integral calculations without the use of any ad- 
justable parameter. Even when we exclude the effects of 
air, viscosity, and surface tension, the collapse is found 
to be not self-similar in a strict sense, but governed by 
power laws with non-universal, Froude-dependent expo- 
nents. Self-similarity is recovered only in the limit of 
infinite Froude number, where the influence of gravity 
becomes negligible and the collapse is truly inertially 
driven. 
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